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Fusion rules in turbulence specify the analytic structure of many-point correlation functions of the 
turbulent field when a group of coordinates coalesce. We show that the existence of flux equilibrium 
in fully developed turbulent systems combined with a direct cascade induces universal fusion rules. 
In certain examples these fusion rules suffice to compute the multiscaling exponents exactly, and in 
other examples they give rise to an infinite number of scaling relations that constrain enormously 
the structure of the allowed theory. 


In a series of recent papers elements of the analytic 
theory of Navier-Stokes turbulence and passive- 

scalar turbulent advection were presented. In 

this Letter we explain that the structure of the essential 
part of these theories is economically summarized by a 
set of “fusion rules” that determine the analytic struc¬ 
ture of n-point correlation functions when a group of co¬ 
ordinates coalesces. We show here that the fusion rules 
can be deduced from very few general assumptions about 
the nature of the universal flux equilibrium that exists in 
fully developed turbulent systems. Of course, the same 
fusion rules can be also established by direct calculations 
in specific examples. We first deduce the fusion rules, 
then we exemplify their utility in determining scaling ex¬ 
ponents, and lastly we demonstrate how in one example 
the fusion rules follow from first principles. 

Consider a turbulent field u(r, t) which is either a vec¬ 
tor or a scalar and denote the difference w(ri|r 2 ,t) = 
u(r 2 , t) — u(ri, t). We discuss the statistical properties of 
the turbulent held in terms of simultaneous many-point 
generalized structure functions of the types 


tion on the characteristic small scale ij. Fully developed 
turbulence is associated with the limit L/r; —> oo. 

We invoke two assumptions of the Kolmogorov type: 
1 . Universality of the hne scale structure of turbulence; 
by this we mean that the correlation functions of the 
type (i. i) have a universal functional dependence on 
their arguments as long as all the separation distances 
involved are much smaller than L. This means that we 
can £x an arbitrary set of velocity differences on the scale 
of L, and the correlation functions will depend on their 
precise choice only via an overall factor that depends on 
the L-scale motions. Mathematically this is expressed as 
the following property of the conditional average: 


(w(ro|ri)w(ro|r 2 )... w(ro|r„) w(ro|Ri)w(ro|R 2 ) (3) 

.. .w(ro|Rjv)) = 5„(ro|ri,r2 ... r„)$„,w(ro|Ri • ■ -Rw) 


for |Ri — ro| ~ L and jr^ — ro| L. 

2. Scale invariance: all the correlation functions are ho¬ 
mogeneous functions of their arguments in the core of the 
inertial interval ij -€1 |ri — ro| L: 


5„(ro|ri,r2 ... r„) = (w(ro|ri) w(ro|r 2 ).. .w(ro|r„)), (1) 
‘5n,m(ro,ro|ri •. .r„;r„+i.. .r„+m) = (w(ro|ri) (2) 

X w(ro|r 2 )... w(ro|r„)w(rQ|r„+i)... w(rQ|r„+m)) , 

etc. Note that when u is a vector the n-point correlation 
is an n-rank tensor. In other words these are correlation 
functions of differences with respect to one, two or more 
reference points. The class of systems that we discuss 
are driven on a characteristic scale referring as the outer 
scale L. This driving can be achieved by either a time 
dependent low frequency “stirring force” or by specify¬ 
ing given values of u at a set of “boundary” points with 
a characteristic separation L away from our observation 
points ro,rg, etc. The systems have dissipation (viscos¬ 
ity, diffusivity etc.) and in the dissipationless limit there 
exists an integral of motion which we refer to as “energy”. 
We consider systems with a “direct” cascade in which the 
intake of energy on the large scales is balanced by dissipa¬ 


5„(Aro|Ari, Ar 2 ... Ar„) = A^'‘5„(ro|ri, r 2 .. .r„), (4) 

where (n is the scaling exponent of the n-order structure 
function. We are particularly interested in systems in 
which is a nonlinear function of n. We refer to such 
systems as “multiscaling”. 

The derivation of these two properties from hrst prin¬ 
ciples differs from system to system. In this Letter we 
discuss the fusion rules and their consequences in sys¬ 
tems for which these assumptions are valid. The hrst set 
of fusion rules that we derive concerns when p points 
(p < n) coalesce with rg, (so that the typical separa¬ 
tion from To is r) and all the other separations remain 
much larger, of the order oi R, r R L. Without 
loss of generality we can choose these p coordinates as 
ri, r 2 ,... Tp. We claim that 

5„(ro|ri,r2,... ,r„) 

= 5p(ro|ri,r2,...,rp)^'„,p(ro|rp+i,rp+2,...,r„) , (5) 


1 



where 4'n^p(ro|rp_|_i, rp_|_ 2 ,..., r„) is a homogeneous func¬ 
tion with a scaling exponent (n — Cp- The derivation of 
the fusion rule (||) follows from Bayes’ theorem and as¬ 
sumptions 1,2. We write 


5„(ro|ri,r2,... ,r„) = j dw(ro|rp+i)... dw(ro|r„) 

X w(ro|rp+i)... w(ro|r„)P[w(ro|rp+i .. .w(ro|r„)] 

X (w(ro|ri), w(ro|r 2 )... w(ro|rp) w(ro|rp+i) 

X w(ro|rp+ 2 )... w(ro|r„)) , (6) 


where 7^[w(ro|rp_|_i ... w(ro|r„)] is the probability to see 
the tensor w(ro|rp+i.. .w(ro|r„). Next comes the cen¬ 
tral idea of this Letter: the properties of flux equilibrium 
and the universal structure of the correlation functions on 
the scale r are the same independen of whether we force 
the system on the scale L ^ r or on the scale R ^ r. 
The conditional average in is proportional to Sp, and 
hence (||). This result seems rather obvious at this point, 
but we will see that it leads to a totally unconventional 
scaling structure of the theory. We should stress that for 
Navier-Stokes and passive scalar advection these fusion 
rules for p = 2 were derived from first principles . 

The next set of fusion rules is obtained for the struc¬ 
ture function Sn,m of when two groups of p < n and 
q < m points coalesce onto Tq and Tq respectively. The 
separation between the groups of point is large and of 
the order of R. The derivation of the fusion rules is now 
obvious, with the result 

^ n,Tn (1*0 7 ^0 1^1 ? ■ ■ ■ 5 ^ 1 * 71+1 j • ■ ■ j 1*77+777) 

= 5p(ro|ri, r 2 ,..., rp)5,(ro|ri, r 2 ,..., r,) 

^ ^1^77,777,2?,^ (1*0 7 ^0 I^P+l ’ ■ ■ ■ ’ ’ ^^+1+1 ’ ■ ■ ■ ’ + * (^) 


The scaling exponent of is C 77+777 - Cp - Cg- Note 

that the fusion rules (||) and (Q) are not decompositions 
into products of lower order correlation functions, and 
the functions 'k are not correlations of velocity differences 
across large separations. In fact we will show that T is 
much larger than the corresponding correlation functions 
in all situations with multiscaling. Evidently one can de¬ 
rive similar fusion rules for three, four or more groups 
of coalescing points with large separations between the 
groups. The structure of the resulting correlation func¬ 
tion will be a product of the correlation function associ¬ 
ated with each group times some function of big sep¬ 
arations which carries the overall exponent. 

Next we discuss fusion rules for correlation functions 
that include gradient fields. These rules depend on the 
type of rotational invariant that one can define from the 
tensors that appear after taking gradients. We will only 
consider the lowest order invariant which is a scalar under 
rotation. For passive scalars T this is |VT • VTp and 
for a vector u the quantity |Vup is the square of the 
strain tensor SijSij where Sij = {dui/drj + duj/dri) j2. 
Consider the quantity 


T2p,77(ro|r2p+i.. .r„) = (|Vu(ro)|^^ 

X w(ro|r 2 p+i)... w(ro|r„)) . (8) 


To evaluate J 2 p,n we consider a related object in which 
all the gradients are taken at different points: 


J2p,n = vil...v:.iviic 


*1 ,Ji- ••'^p ^3p 

fei ,Zi.. .fep ,jp 


X (ty'=i(ro|ri)'u;'i(ro|ri)... ■u;'’’p(ro|rp)w''>(ro|rp 
X w(ro|r 2 p+i)... w(ro|r„)). 


(9) 


where the contraction tensor C ensures that the 
required scalar is obtained. We will represent 
this quantity in a compact form without display¬ 
ing all the tensor indices as J 2 p, 77 (ro|ri, r 2 ... r„) = 
VriVrj ... Vr^C5„(ro|ri,r'^.. .r„). The quantity (|) 
gives us J 2 p ,77 when the 2p first points coalesce together 
with To, whereas all the rest of the coordinates remain 
a typical distance R from Tq. When R is in the inertial 
interval we expect scaling behaviour in terms of R, 

J2p,n OC _ (10) 


Considering the distances between all the coalescing 
points to be in the inertial range we apply and find 
for 2p coalescing points 

T2p,77(ro|i*i,r2 .. .r„) = V,-! ... Vr^52p(ro|ri, r'^ .. .r),) 

X ^'77,2p(i*o|r2p+i,r2p+2,... ,r„) . (11) 


We expect that J 2 p,n is independent of the first 2p coor¬ 
dinates when the distances between them are well in the 
viscous regime. Our next fundamental assumption is that 
there exists a characteristic viscous length r]{2p,n, R) 
at which J 2 p,n crosses smoothly from inertial range be¬ 
haviour to dissipative behaviour with respect to the first 
2p coordinates. This allows us to evaluate the coalescing 
gradients by taking the 2p separations to be ri(2p, n, R): 

T2p,77(ro|ri,r2,... ,r„) 'q{2p,n, R)‘^'^p~'^p (12) 

X ^' 77 , 2 p(i*o|i* 2 p+i, i* 2 p+ 2 , • ■ •, Tri) (2p coalescing points). 

If there are two groups of coalescing points with gradi¬ 
ents, with p points coalescing onto rg and q points coa¬ 
lescing on Tq respectively, we consider Jp^q^n.m (where as 
before n + m>p + q is the total number of points). The 
rule for p and q coalescing points is 

Tp,9,77,777(i*o|i*i, r 2 ,..., r„) ~ ? 7 (p, n, R)^^~‘Pr]{q, n, 

^ ^n,m,p,q (ro|rp+i,rp+ 2 ,...,r„) . (13) 

The generalization of this fusion rule for three or more 
groups of coalescing points with gradients is obvious. 

This is as much as one can do in general. Now the 
crucial issue is how rj{2p, n, R) depends on its arguments. 
The simplest version of the theory comes about when 
the dissipative length is independent of R, r]{2p, n, R) = 
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r]{2p,n). This is realized for example in passive scalar 
convection as is shown below. In this case the fusion 
rules imply various sets of scaling relations. For example 
the exponents ^{2p,n) of J 2 p,n are given by 

C(2p,n) = Cri - C2p • (14) 

As another example of scaling relations consider the cor¬ 
relation functions 

= <;|Vu(r)|2«|Vu(r + R)l2*) oc . (15) 

From (|^) in the case n = m = p = q = 2s we get a set 
of scaling relations 

p{2s) = 2 ( 2 s - C 4 s • (16) 


Next we can consider a correlation of I gradient fields 
with the same power, (i.e. |Vup®) at I different points 
which are separated by a distance of the order of R. The 
corresponding exponent p{l, 2s) is 

p{l,2s)=K2s-(:2sl ■ (17) 

This algebra can be generalized to any correla¬ 
tion of powers of |Vu|^. For example, the scal¬ 
ing exponent p,{pi,p 2 , ■ ■ -Pn) of a correlation of fields 
(|Vu(ri)|Pi |Vu(r 2 )|^^ ... |Vu(r„)|P") in which all the 
separations is of the order of R is 


p{pi,P2,-■-Pn) ='^Cpj - Cp , P = '^Pj ■ (18) 

f=i 

In usual operator algebras [p^ , pd]|l^Jl^ every local field 
is associated with a leading exponent and the correla¬ 
tion function scales with the sum of these exponents. In 
this case the algebra is different. There is a global expo¬ 
nent (^p from which one subtracts the sum of individual 
exponents Cpj ■ In multiscaling situations the global ex¬ 
ponent is a nonlinear function of p = '^Pj- Accordingly 
it is not a property of every individual field. We note 
here without demonstration that invariants of the gradi¬ 
ent field tensors other than scalars are associated with 
other individual exponents. 

The range of applicability of these fusion rules should 
be understood on the basis of the equations of motion 
for any given system. As an example we explain here 
briefly why they are applicable for Kraichnan’s model 0 
of passive scalar convection with a driving velocity field 
that is (5-correlated in time. It was shown that the 
cumulant part of the 2 n-order correlation function 
F 2 n = {T{vi,t)T{Y 2 ,t).. .T{v 2 n,t)) satisfics for n > I 
the exact homogeneous differential equation 


2 n 


a 




2n 


7 ^ 2 «(ri,i‘ 2 ,...,r 2 n) = 0 , (19) 


where k is the molecular diffusivity and is the Lapla- 
cian operator acting on Tq, . The operator B 2 n is the sum 
of the binary operators Bap: 


BafS — ^/S) 






B2n = 


2n 

E 

a>/3=l 


( 20 ) 


Here /iij (R) is the eddy diffusivity that behaves like 
Hi?*’'* with 0 < ()h < 2 .The scaling exponent C 2 satis¬ 
fies 1 ^ the exact relation (,2 = 2 — (h- 

In the inertial range of scales we can disregard the 
Laplacian operators in this equation. For deriving the 
fusion rules (|^) we consider the p coalescing points with 
characteristic separation r and denote their coordinates 
by the index a oi a'. The remaining 2n — p coordinates 
have characteristic separations R and are denoted by (3 or 
B'. We assemble 0 the B operators in three groups: Bp = 

Ea>a' ^2n-p = E^>/3' and B^ = Y^ap^^P- 

The evaluation of the action of the operators in the first 
and second groups is H/r‘'^ and H/R^^ respectively. The 
evaluation of the action of each term in the third group 
is HR/rR'^^. However space homogeneity results in a 
cancellation of this evaluation in the sum of the terms in 
this group. The next order surviving evaluation is again 
H/R'^^. We thus combine the second and third group 
into an effective operator B. The equation to consider is 

[Bp + B]F2n = 0. (21) 

When p = 2 we can find the solution of ( |^ ) as the fol¬ 
lowing expansion in powers of the small difference ri 2 : 
A 2 {R} + r‘llC 2 {R} + rl^^D 2 {R} + rl 2 E 2 {R} +■ ■ •, where 
A 2 , C 2 , D 2 etc. are some functions of the large separa¬ 
tions of the order of R. When we use this solution to 
compute S 2 n the leading contribution A{i?} drops and 
we find (||) for p = 2. For p > 2 we need to distinguish 
between even and odd p’s because of the special property 
of passive advection in which iS 2 n-i-i = 0. The next even 
p is p = 4. For this case we find a solution in the form 


4 

FY = A,{R} + C,{R}[ ^ 

OLOi' — l 


( 22 ) 


-I- D4{R}[{ri2ri3)'^'^ + {ri2ri4)^^ + (ri2r23)^’' -f .. . ] 
-I-F'Kri, r 2 , ra, r4)'I'2n,4{7?} + ... , 


where F 4 is a contribution of a new type, as it solves the 
homogeneous equation ( p^ . Computing S 2 n the first two 
terms disappear and in a multiscaling situation the lead¬ 
ing contribution becomes the last. In fact, this is the 
general rule for any even order, and is the explicit mech¬ 
anism for the fusion rules in this particular case. It arises 
here from the possibility to split the total operator B 2 n 
into the two groups Bp and B such that for p coalescing 
points Bp carries the leading contribution. Since the sum 
of Laplacians in ( 0 ) is also dominated by the sum up 
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to p, the crossover scale rj{p, 2n, R) from inertial range 
to dissipative behaviour is determined in this case by a 
balance between — Bp. It therefore can¬ 

not depend on n or on R\ ri{p, 2n, R) = r]{p). The fusion 
rules (p^-(p^ which were based on the independence of 
r] on R follow. 

In fact, these results, and in particular the scaling re¬ 
lations ( p^ seem sufficient to determine the exponents 
in their entirety. The necessary steps were detailed in 
and will not be repeated here. 

The case of Navier-Stokes turbulence calls for addi¬ 
tional considerations. The fusion rules (||), (|^) were found 
from first principles for p = 2 Q and we believe that 
similar techniques can be used to establish them for any 
p. Eqs.(|l^- (^) follow, but in the Navier-Stokes case 
it is possible that the dissipative scale r]{p, n, R) is R 
dependent. If we assume that this is not the case the 
scaling relations obtained above should apply also to the 
Navier-Stokes case. The consequences of such an assump¬ 
tion were discussed in detail in j|]. To explore another 
possibility we follow Kolmogorov’s refined similarity hy¬ 
pothesis [p^ in assuming that the conditional average 
1^(1 Vup|w(0|R)) ~ u>(0|R)^/i?. This assumption means 


uJ2MR}--Sn+i{R)/R. (23) 

Comparing with (H this can be consistent only if 


■77(2, n, i?)! 2-C2 / R\Cn-Cn-l+C 3 -C 2 

. 77(2) J VlJ 


(24) 


where 77 ( 2 ) is by definition the viscous cutoff for the sec¬ 
ond order structure function, 7252 ( 77 ( 2 )) ^ r]{2)‘^e. The 
Holder inequalities guarantee that Cn — Cn-i is a decreas¬ 
ing function of n in a multiscaling situation. Accordingly, 
the effective dissipative scale of J 2 ,n for two coalescing 
points 77 ( 2 , 77 , ii) is much smaller than the viscous cutoff 
for 52 , 77 ( 2 ). 

Needles to say, with this assumption all our scaling 
relations change. For example consider ^(2, n) of (10). 
Instead of ([IJ) we have now 


e( 2 , n) = Cn+i - 1 . 
Another example is KlJ. we find now 


(25) 


p(l) = 2 - Ce . (26) 


This result is known as “the bridge” in the phenomeno¬ 
logical theory of multiscaling turbulence, c.f. |p^ , p^ . 
Notwithstanding the different scaling relation, the oper¬ 
ator algebra that is induced has “global” and individual 
scaling exponents as discussed above. The values of these 
exponents may be changed due to the R dependence of 
the dissipative cutoff as it appears in different models. 

In summary, we proposed fusion rules for multiscaling 
turbulent statistics that induce an unusual operator al¬ 
gebra. These rules are of two classes. The first does not 


involve gradients and is universal for all turbulent sys¬ 
tems with a direct cascade of “energy” in which there 
exists a universal flux equilibrium. The second class in¬ 
volves gradients and these bring in an explicit depen¬ 
dence on a viscous scale that in general is not universal. 
We explained why in the case of passive scalar advection 
this problem may be solved. Accordingly one can derive 
an infinite set of non trivial scaling relations that allow 
the expression of the scaling exponents of the correla¬ 
tion functions of gradient fields with the exponents (^n 
of the structure functions. For Navier-Stokes turbulence 
the fusion rules that involve gradients must by supple¬ 
mented with a theory for the R dependence of the vis¬ 
cous cutoff. We exemplified the influence of a reasonable 
assumption about this dependence, but a solid theory of 
this dependence based on the Navier-Stokes equations is 
still a future project. 
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